
HL Paper 3
This question explores models for the height of water in a cylindrical container as water drains out.
The diagram shows a cylindrical water container of height 3.2 metres and base radius 1 metre. At the base of the container is a small circular valve, which enables water to drain out.
Eva closes the valve and fills the container with water.
At time t=0, Eva opens the valve. She records the height, h metres, of water remaining in the container every 5 minutes.
Eva first tries to model the height using a linear function, h(t)=at+b, where a, b∈ℝ.
Eva uses the equation of the regression line of h on t, to predict the time it will take for all the water to drain out of the container.
Eva thinks she can improve her model by using a quadratic function, h(t)=pt2+qt+r, where p, q, r∈ℝ.
Eva uses this equation to predict the time it will take for all the water to drain out of the container and obtains an answer of k minutes.
Let V be the volume, in cubic metres, of water in the container at time t minutes.
Let R be the radius, in metres, of the circular valve.
Eva does some research and discovers a formula for the rate of change of V.
dVdt=-πR2√70 560h
Eva measures the radius of the valve to be 0.023 metres. Let T be the time, in minutes, it takes for all the water to drain out of the container.
Eva wants to use the container as a timer. She adjusts the initial height of water in the container so that all the water will drain out of the container in 15 minutes.
Eva has another water container that is identical to the first one. She places one water container above the other one, so that all the water from the highest container will drain into the lowest container. Eva completely fills the highest container, but only fills the lowest container to a height of 1 metre, as shown in the diagram.
At time t=0 Eva opens both valves. Let H be the height of water, in metres, in the lowest container at time t.
Find the equation of the regression line of h on t.
Interpret the meaning of parameter a in the context of the model.
Suggest why Eva’s use of the linear regression equation in this way could be unreliable.
Find the equation of the least squares quadratic regression curve.
Find the value of k.
Hence, write down a suitable domain for Eva’s function h(t)=pt2+qt+r.
Show that dhdt=-R2√70 560h.
By solving the differential equation dhdt=-R2√70 560h, show that the general solution is given by h=17 640(c-R2t)2, where c∈ℝ.
Use the general solution from part (d) and the initial condition h(0)=3.2 to predict the value of T.
Find this new height.
Show that dHdt≈0.2514-0.009873t-0.1405√H, where 0≤t≤T.
Use Euler’s method with a step length of 0.5 minutes to estimate the maximum value of H.
Alessia is an ecologist working for Mediterranean fishing authorities. She is interested in whether the mackerel population density is likely to fall below 5000 mackerel per km3, as this is the minimum value required for sustainable fishing. She believes that the primary factor affecting the mackerel population is the interaction of mackerel with sharks, their main predator.
The population densities of mackerel (M thousands per km3) and sharks (S per km3) in the Mediterranean Sea are modelled by the coupled differential equations:
dMdt=αM-βMS
dSdt=γMS-δS
where t is measured in years, and α, β, γ and δ are parameters.
This model assumes that no other factors affect the mackerel or shark population densities.
The term αM models the population growth rate of the mackerel in the absence of sharks.
The term βMS models the death rate of the mackerel due to being eaten by sharks.
Suggest similar interpretations for the following terms.
An equilibrium point is a set of values of M and S , such that dMdt=0 and dSdt=0.
Given that both species are present at the equilibrium point,
The equilibrium point found in part (b) gives the average values of M and S over time.
Use the model to predict how the following events would affect the average value of M. Justify your answers.
To estimate the value of α, Alessia considers a situation where there are no sharks and the initial mackerel population density is M0.
Based on additional observations, it is believed that
α=0.549,
β=0.236,
γ=0.244,
δ=1.39.
Alessia decides to use Euler’s method to estimate future mackerel and shark population densities. The initial population densities are estimated to be M0=5.7 and S0=2. She uses a step length of 0.1 years.
Alessia will use her model to estimate whether the mackerel population density is likely to fall below the minimum value required for sustainable fishing, 5000 per km3, during the first nine years.
γMS
δS
show that, at the equilibrium point, the value of the mackerel population density is δγ;
find the value of the shark population density at the equilibrium point.
Toxic sewage is added to the Mediterranean Sea. Alessia claims this reduces the shark population growth rate and hence the value of γ is halved. No other parameter changes.
Global warming increases the temperature of the Mediterranean Sea. Alessia claims that this promotes the mackerel population growth rate and hence the value of α is doubled. No other parameter changes.
Write down the differential equation for M that models this situation.
Show that the expression for the mackerel population density after t years is M=M0eαt
Alessia estimates that the mackerel population density increases by a factor of three every two years. Show that α=0.549 to three significant figures.
Write down expressions for Mn+1 and Sn+1 in terms of Mn and Sn.
Use Euler’s method to find an estimate for the mackerel population density after one year.
Use Euler’s method to sketch the trajectory of the phase portrait, for 4≤M≤7 and 1.5≤S≤3, over the first nine years.
Using your phase portrait, or otherwise, determine whether the mackerel population density would be sufficient to support sustainable fishing during the first nine years.
State two reasons why Alessia’s conclusion, found in part (f)(ii), might not be valid.
This question explores methods to determine the area bounded by an unknown curve.
The curve y=f(x) is shown in the graph, for 0⩽.
The curve passes through the following points.
It is required to find the area bounded by the curve, the -axis, the -axis and the line .
One possible model for the curve is a cubic function.
A second possible model for the curve is an exponential function, , where .
Use the trapezoidal rule to find an estimate for the area.
With reference to the shape of the graph, explain whether your answer to part (a)(i) will be an over-estimate or an underestimate of the area.
Use all the coordinates in the table to find the equation of the least squares cubic regression curve.
Write down the coefficient of determination.
Write down an expression for the area enclosed by the cubic function, the -axis, the -axis and the line .
Find the value of this area.
Show that .
Hence explain how a straight line graph could be drawn using the coordinates in the table.
By finding the equation of a suitable regression line, show that and .
Hence find the area enclosed by the exponential function, the -axis, the -axis and the line .
This question will investigate the solution to a coupled system of differential equations and how to transform it to a system that can be solved by the eigenvector method.
It is desired to solve the coupled system of differential equations
where and represent the population of two types of symbiotic coral and is time measured in decades.
Find the equilibrium point for this system.
If initially and use Euler’s method with an time increment of 0.1 to find an approximation for the values of and when .
Extend this method to conjecture the limit of the ratio as .
Show how using the substitution transforms the system of differential equations into .
Solve this system of equations by the eigenvalue method and hence find the general solution for of the original system.
Find the particular solution to the original system, given the initial conditions of part (b).
Hence find the exact values of and when , giving the answers to 4 significant figures.
Use part (f) to find limit of the ratio as .
With the initial conditions as given in part (b) state if the equilibrium point is stable or unstable.
If instead the initial conditions were given as and , find the particular solution for of the original system, in this case.
With the initial conditions as given in part (j), determine if the equilibrium point is stable or unstable.
This question will investigate the solution to a coupled system of differential equations when there is only one eigenvalue.
It is desired to solve the coupled system of differential equations
The general solution to the coupled system of differential equations is hence given by
As the trajectory approaches an asymptote.
Show that the matrix has (sadly) only one eigenvalue. Find this eigenvalue and an associated eigenvector.
Hence, verify that is a solution to the above system.
Verify that is also a solution.
If initially at find the particular solution.
Find the values of and when .
Find the equation of this asymptote.
State the direction of the trajectory, including the quadrant it is in as it approaches this asymptote.
This question is about modelling the spread of a computer virus to predict the number of computers in a city which will be infected by the virus.
A systems analyst defines the following variables in a model:
- is the number of days since the first computer was infected by the virus.
- is the total number of computers that have been infected up to and including day .
The following data were collected:
A model for the early stage of the spread of the computer virus suggests that
where is the total number of computers in a city and is a measure of how easily the virus is spreading between computers. Both and are assumed to be constant.
The data above are taken from city X which is estimated to have million computers.
The analyst looks at data for another city, Y. These data indicate a value of .
An estimate for , can be found by using the formula:
.
The following table shows estimates of for city X at different values of .
An improved model for , which is valid for large values of , is the logistic differential equation
where and are constants.
Based on this differential equation, the graph of against is predicted to be a straight line.
Find the equation of the regression line of on .
Write down the value of , Pearson’s product-moment correlation coefficient.
Explain why it would not be appropriate to conduct a hypothesis test on the value of found in (a)(ii).
Find the general solution of the differential equation .
Using the data in the table write down the equation for an appropriate non-linear regression model.
Write down the value of for this model.
Hence comment on the suitability of the model from (b)(ii) in comparison with the linear model found in part (a).
By considering large values of write down one criticism of the model found in (b)(ii).
Use your answer from part (b)(ii) to estimate the time taken for the number of infected computers to double.
Find in which city, X or Y, the computer virus is spreading more easily. Justify your answer using your results from part (b).
Determine the value of and of . Give your answers correct to one decimal place.
Use linear regression to estimate the value of and of .
The solution to the differential equation is given by
where is a constant.
Using your answer to part (f)(i), estimate the percentage of computers in city X that are expected to have been infected by the virus over a long period of time.
Consider the system of paired differential equations
.
This system is going to be solved by using the eigenvalue method.
If the system has a pair of purely imaginary eigenvalues
Show that if the system has two distinct real eigenvalues then .
Find two conditions that must be satisfied by , , , .
Explain why and must have opposite signs.
In the case when there is a pair of purely imaginary eigenvalues you are told that the solution will form an ellipse. You are also told that the initial conditions are such that the ellipse is large enough that it will cross both the positive and negative axes and the positive and negative axes.
By considering the differential equations at these four crossing point investigate if the trajectory is in a clockwise or anticlockwise direction round the ellipse. Give your decision in terms of and . Using part (b) (ii) show that your conclusions are consistent.
Find the value of .
Illustrate graphically the inequality .
Hence write down a lower bound for .
Find an upper bound for .
The number of brown squirrels, , in an area of woodland can be modelled by the following differential equation.
, where
One year conservationists notice that some black squirrels are moving into the woodland. The two species of squirrel are in competition for the same food supplies. Let be the number of black squirrels in the woodland.
Conservationists wish to predict the likely future populations of the two species of squirrels. Research from other areas indicates that when the two populations come into contact the growth can be modelled by the following differential equations, in which is measured in tens of years.
, , ≥ 0
, , ≥ 0
An equilibrium point for the populations occurs when both and .
When the two populations are small the model can be reduced to the linear system
.
For larger populations, the conservationists decide to use Euler’s method to find the long‑term outcomes for the populations. They will use Euler’s method with a step length of 2 years ().
Find the equilibrium population of brown squirrels suggested by this model.
Explain why the population of squirrels is increasing for values of less than this value.
Verify that , is an equilibrium point.
Find the other three equilibrium points.
By using separation of variables, show that the general solution of is .
Write down the general solution of .
If both populations contain 10 squirrels at use the solutions to parts (c) (i) and (ii) to estimate the number of black and brown squirrels when . Give your answers to the nearest whole numbers.
Write down the expressions for and that the conservationists will use.
Given that the initial populations are , , find the populations of each species of squirrel when .
Use further iterations of Euler’s method to find the long-term population for each species of squirrel from these initial values.
Use the same method to find the long-term populations of squirrels when the initial populations are , .
Use Euler’s method with step length 0.2 to sketch, on the same axes, the approximate trajectories for the populations with the following initial populations.
(i) ,
(ii) ,
Given that the equilibrium point at (800, 600) is a saddle point, sketch the phase portrait for ≥ 0 , ≥ 0 on the same axes used in part (e).